Derivation of Hyperbolic Transfer Equations from BGK-Equation 
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We use the integral form of the Boltzmann equation which aUows us to take into account the 
memory effects using the initial condition that selects the solutions going to the local equilibrium 
Maxwell distribution in the t — > — oo limit. Implementing the relaxation-time approximation for 
the collision integral (BGK-equation) we present the derivation of the hyperbolic Navier-Stokes and 
the hyperbolic heat conduction equations in the first order approximation. It is shown that the 
relaxation time in the obtained hyperbolic equations is the Maxwellian relaxation time. As special 
case we obtain the telegraph equation for the heat propagation in static medium and estimate the 
relaxation time for the heat conduction in some materials. 

PACS numbers: 



I. INTRODUCTION 

The hyperbolic equation for the heat conduction (and 
also for the diffusion) has been of interest for a long time. 
It is known that this equation has so-called the relaxation 
term which takes into account the propagation of a signal 
with a finite speed. Many authors have dealt with this 
problem. Cattaneo 0, Q and Vernotte Q were the first 
who based their analysis on kinetic theory arguments and 
added the relaxation term to the Fourier equation. The 
Cattaneo equation or the hyperbolic heat flux has the 
form 



(1) 



where A is the thermal conductivity and Tg is some re- 
laxation time. In our designations Tg is the hyperbolic 
relaxation time, which is usually obtained from the ex- 



periment. For example, for nitrogen Tg 



10- 



and 



for aluminium Tg = 10~^^ s |^. From Eq. Q and the 
continuity equation one can write the hyperbolic heat 
conduction equation. 

This point of view has been used in several contexts: 
e.g. in the analysis of waves in thermoelastic media, fast 
explosions, and a second sound in solids. Wide ranging 
bibliography can be found in 0,1113,1111. An excellent 
monograph on this problem is in Ref. [lOj. 

The problem of derivation of the hyperbolic equation 
from microscopic equations caused the special interest 
for a long time. There are several works on the hyper- 
bolic diffusion. One of the first publications appeared in 
1935 pT|. In this paper Davidov obtained the hyperbolic 
equation only for one-dimensional case and for a constant 
molecular velocity. The similar results on this problem 
one can find in other works [T^ IT^ [l^ ITsl IT^ . In the 
publications of Davies |0| and Das |l8l| on Brownian mo- 
tion the hyperbolic diffusion equation was obtained from 
the Fokker-Plank equation. It means that equations for 
transfer processes of the hyperbolic form can be derived 
from the underlying kinetic equations. 



For the heat transfer the Grad's model is of great in- 
terest and its thirteen- moment approximation. This 
theory solves the Boltzmann equation in non-equilibrium 
situations, and the hyperbolic heat flux and others 
hyperbolic-like transfer equations can be obtained in the 
linear approximation [20|. It is shown that the relax- 
ation time in these equations explicitly depend on the 
collision parameters. However, exact form of Eq. Q has 
not been obtained using Grad's method. Instead of Eq. 
([T|) the hyperbolic heat flux in Grad's expansion contains 
other terms that can be explained by the expansion of the 
collision integral. 

Nevertheless, there is another way for the theoretical 
investigation of the non-equilibrium processes by con- 
structing the distribution function / which dep ends on 
time in the form of a delay functional 0, l22l 123^. i.e. 
taking into account memory effects. There is a reduced 
description of the non-equilibrium system where the dis- 
tribution function / is expressed in terms of the first 
five moments, which are parameters of the hydrodynamic 
state (concentration n, mass velocity u and temperature 
T). This method is used in the well-known Chapman- 
Enskog approach 0; EE Ell- The function / is sought as 
a functional of the latter parameters which satisfy the hy- 
drodynamic equations. These solutions are known as nor- 
mal solutions [l^ I23L |2J| . In contrast to the Chapman- 
Enskog approach, in the work of Zubarev the Boltz- 
mann equation is expressed in an integral form, where 
there is the initial condition that provides the selection 
of the normal solutions. Consequently it has no other 
solutions. The analogous method was used for the Liou- 
ville equation in several works 23, 27, 28, 29, 30] . For 
the derivation of the transfer equations with physically 
realistic time delays one more condition is necessary such 
as, for example, the relaxation time approach. If we do 
not use any additional supposition, it is very difficult to 
choose the physical relaxation time explicitly, while Eq. 
([T|l is approximation, and as it will be shown in this work, 
the relaxation time approximation. 

In the present paper we obtain the hyperbolic Navier- 
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Stokes and heat conduction equations from the Boltz- 
mann equation in the relaxation time approximation 
(BGK-equation j^) using the idea of the delay func- 
tional. This approach allows to obtain the expressions 
for the relaxational fluxes in integral as well as in dif- 
ferential forms. For this result we use the integral form 
of the Boltzmann equation with an initial condition, and 
the relaxation time approximation (RTA) for the collision 
integral. 

In Sec. ^we give a short description of the basic prop- 
erties of the Boltzmann equation. In Sec. IIIII we con- 
sider the integral form for the Boltzmann equation and 
the relaxation time approximation for the collision inte- 
gral. The expression for the pressure tensor, heat flux 
and hyperbolic balance equations are obtained in Sec. 
IIVI Then we consider the hyperbolic heat propagation 
in static medium in Sec. Conclusions are given in Sec. 

ED 

For the illustration we consider mono-atomic gas and 
distribution function in the first approximation with re- 
spect to small gradients. 



where n{r,t) is the particle number density, u{r,t) is 
the mass velocity (or hydrodynamic velocity), T{r,t) is 
the temperature, c — v u{r,t) is the thermal velocity, 
k is the Boltzmann constant, and m is the mass of the 
particle. 

Multiplying Eq. ||2J) by 1, Va, and and integrating 
with respect to v we obtain the known balance equations 



dn d 



+ -7^ — {nua) = 0, 



dt dva 



dUa 



up 



I dp, 



dUa 

drp p drp ' 



'q/3 



(5a) 



(5b) 



where p — mn is the mass density. The pressure tensor 
Pap and the heat flux are defined by the relations 



II. THE BASIC PROPERTIES OF THE 
BOLTZMANN EQUATION 



Pa/3 = / mCaCi3f{r,V,t)(fv, 



(6a) 



Let us consider the Boltzmann equation for the single- 
particle distribution function of the molecules over veloc- 
ities V in a, point r at an instant time t 



Ca—^j[r,v,t)(Fv. 



(6b) 



(2) 



III. INTEGRAL FORM FOR THE BOLTZMANN 
EQUATION 



The collision integral is given by 

/(./,/) = I {f[f ^ fif) 9 <T d'n d'v,, 

f[ = f{r,v[,t), f = f{r,v\t), 
h = .fir,vi,t), f = f{r,v,t), 
g=\v- vi\ 



(3) 



where v, Vi is the initial velocities, v' , v[ is the final 
velocities of molecules after pair collision, cr is the scat- 
tering cross section, and is the solid angle . 

The parameters of the hydrodynamic state of the 
system are expressed through the distribution function 
f{r,v,t) as 



n{r,t) - / f{r,v,t)d\, 



(4a) 



n{r,t) 



Vaf{r,v,t)d^v, (4b) 
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Since the Boltzmann equation is very complicated, 
simpler kinetic equations have been proposed in the liter- 
ature. A widely used simplified model is the well known 
the relaxation time approximation (RTA) for the colli- 
sion integral (e. g. see Refs. [Toll25ll2^l3lll3^ '). In this 
case the resulting Boltzmann equation has the form 



dj_ 

dt 



where 



(0) 



1 2 p 2kT ^ 



(7) 



(8) 



is the local equilibrium Maxwell distribution, is the 
relaxation time which is of the order of the mean free 
time tq. 

Let us write the function / in the following form: 



/ = F + /("). 



(9) 



Substituting @ in Eq. ^ we obtain the equation for 
the function F 



dF 1 

dt Tr 



F 



Df(o) dF 



(10) 
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where 



D 
Dt 



d_ 

dt 



is the time derivative along the trajectory of a particle. 

We note that that the hydrodynamic variables n{r,t), 
u{r,t) and T{r,t), defined by Eqs. coincide with 

the corresponding variables calculated from the local 
Maxwellian distribution (jSJ. Thus, the first five moments 
of the function F = f — /^"^ equal zero. 

To specify a solution of Eq. UlUfl the idea of Zubarev 
and Honkin [23| is used in the following considerations. 
One assumes that at some time tg we have F(io) = 0, 
i. e. /(to) — /'"'(^o)- Then there considers the limiting 
process t — to ^ oo supposing that the solution goes to 
the local Maxwell distribution in the t — > — oo limit. To 
satisfy the latter conditions Eq. l(Tn|) is represented in 
the integral form 



Fit) 



Dt '"9rJ„ 



dt', (11) 



where the subscript t' behind the bracket means that the 
time argument inside the bracket equal t' . Thus, Eq. 
(|ll|l allows us to take into account the memory effects, 
and to select normal solutions corresponding to an initial 
condition 



lim F{t) ^ 0. 



(12) 



In an assumption of small gradients of hydrodynami- 
cal parameters the integral equation (|11|) can be solved 
by expanding F in series in the powers of the Knudsen 
number Kn (Kn = Ij/Al, If is the mean free path of 
particle, and Al is the characteristic space scale for the 
hydrodynamic processes), or iteratively, assuming 



(13) 

This itcrational procedure can be implemented for the 
consequent determining of different approximations in 
constructing of normal solutions. 

In fact, one can obtain the same result if we use for- 
mally the Zubarev's method ^ in the RTA. The differ- 
ence is only in the exponential term where the term 1 /ti- 
is substituted by (1 -I- eTr)/Tr. The parameter e one can 
consider as a some parameter providing the convergence 
of Eq. Hll|l . In our case the role of this parameter plays 
the exponential term I/t^ with e = 0. 



IV. THE FIRST APPROXIMATION 

As the first approximation in Eq. H13|l we set F'"^ = 
(i. e. from Eq. © F^°^ = /^°^ - /^°^ = 0). Thus, 



where the gradients of F^^^ are neglected . On the other 
hand by differentiating Eq. (|14|l over time we obtain 



dt Tr Dt ■ 



(15) 



i. e., the first time derivative d/dt of F^^^ is taken into 
account. It can be explained by the fact that the integral 
equation Hll|l does not have the derivative dF/dt, which 
is not considered in the iterational expansion. 



The Transfer Lows 



Using Eq. © we have 



Dt 



RT 



RT 



1 dn dua 
n dt dra 
dua 1 dp 
dt p dva 
2 3\ fdRT 
2, 



2RT 



V dt 



+ 



3 ar„ 



1 2c \dua 



2RT 



where 



d_ 
di 



d_ 

dt 



dRT 
dva 



9r„ 



(16) 



R = k/m is the gas constant. 

As mentioned above, the function F does not give 
contribution to the hydrodynamic parameters. Conse- 
quently, for the first five moments of the function F^^^^ 
the following conditions have to be fulfilled 



F'-^^d^v = 0, 



-F^'^U^v = 0. 



(17a) 



(17b) 



(17c) 



From these conditions there follows that the first three 
terms of Eq. H16|l in round brackets, which are the hy- 
drodynamical equations of ideal fiuid, have to be equal 
zero. To prove this we will follow the work of Zubarev 
and Honkin |2lj | representing their results for our case. 
First, we will show that the integrals 



(18) 



4 



V'/3 




/? = 
/3=1,2,3 



equal zero for a 7^ and do not equal zero for a = (3. 
Let us consider a function 



The function satisfies an equation 



(19) 



(20) 



We multiply Eq. H2U|) on ipa and integrate it over v 
obtaining 



where 



9 1 

^ + — ) ^0/3(0 = "^Q/S, 



(21) 



(22) 



Taking into account that the integrals Jap = for a 7^ /3, 
and Jaf3 7^ for a = /3, we prove the our statement for 
the integrals ((TH|l . 

In similar way one can show that 



d / dt'e^^''-'\f^°'>^ai3)t' =0, (23a) 



d^J dt'e^(*'-*)(/(0)^„),, =0, (23b) 



where 



/ di'e^(*'-*)(/(°Va)t' =0, (23c) 



(24a) 



2RT 



(24b) 



The relations 123|) follow from the properties of the mu- 
tual orthogonality of the functions ipp^ fap and ipa. 

Let us substitute expressions H14|l and H16(l into Eqs. 
(|17|l . Then using the properties of the integrals H18() and 
(|23|l we obtain that the relations (|17|) are fulfilled only 
in the case when the hydrodynamical equations of ideal 
fluid are implemented 



dn 
'dt 



I— 

dra 



(25a) 



dT 
'dt ~ 



1 dp 
pdra 



3 dra 



(25b) 
(25c) 



and the function F^^^ has the form 



-(t'~t) 



RT^'^'d^ 



d RT{t') 
dra 



dt'. (26) 



We see that the right hand side of Eq. is expressed in 
terms of gradients of the hydrodynamic variables. Thus, 
the function F'^^^ represents a normal solution of the 
Boltzmann equation. 

We can obtain the pressure tensor and the heat flux 
in the first approximation using expressions © and or- 
thogonal properties for ipap and Lpa'. 



Pap - P5ap 
dUa _^ 

drp dra 3 dr^ 



M{t, t') X 



dua dup 2 du^ - , ,,, 
X I — -t- ^ - -^5ap \ dt'. 



dra 



(27a) 



(27b) 



where 



L{t,t') 



mc 



d^c 



(28b) 



are the kernels of the equations (|27|l . 

The integral form for the heat flux, like Eq. Ij27b|l . 
is known relatively long time .''>.''>. 34j. These equations 
describe the transfer in the local non-equilibrium systems 
or in the so-called materials with memory [s^. However, 
as shown above the same result one can also obtain for 
the pressure tensor Eq. H27a|) . 

Eqs. (|27|) take the classical form for t 3> r,.. In 
this case the exponent becomes the delta-function, con- 
sequently one can carry out behind the integral the other 
integrands in the point t' = t 



(PaP ~ pSap ] 



^M{t) X 



' dua dup 2 du-y ^ 



drp dra 

, dua dup 2duj \ 
-^^■^'^'\d^+d^-3d^^-^) 



dua du 
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2du^ 



drp dra 3 dr^ 



5ap , 



(29a) 
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, , , dT(t) , dT{t) 
= -TrL{t)—^ = -A— (29b) 



dL{t,t') 
Ft 



m 

YJ [ ' dt 2 



dua (? du 



dt 



CaCf} X 



(34b) 



where fi and A are the classical shear viscosity and ther- 
mal conductivity. Using the Appendix ^ one can show 
that they have forms: 



are the kernels of Eqs. (|33|) . Since the expressions 1)34(1 
contain the orthogonal functions in the meaning of the 
scalar product H18|l . the integrals ((33|1 equal zero 



H — TrM(t) = TrUkT = —Trv'^nm, 



A = TrLit] — Tr TIK = —TrV UK, 

^ ' 2 m 4 



(30a) 



(30b) 



where v = ^j2kT/ m is the average molecular velocity. 
These equations have the known classical limit: 



TrV — > const, if Tr ^ and v oo. 



(31) 



It is important to note that in contrast to the known 
Navier-Stokes and Fourier lows the equations H27() con- 
tain time delay. 



Pal3 =0, qa= 0. 



(35) 



For further considerations we write Eqs. ((32(1 in a more 
convenient form: 



Pal3 



pSap - 2^[ Aq/3 - i:^-^Sa 



where 



3 dr^ 



dt 



^ dT 

dra 



■ dt 



(36a) 



(36b) 



B. Balance Equations 



1 / dua , dup \ 



(37) 



In order to obtain the hydrodynamic balance equations 
in the first approximation we differentiate Eqs. ((27() over 
time: 



dt 



dqg 

dt 



-M{t) 
1 



/ dua dup 2 du-y ^ \ 
\ drp dra 3 dr^ J 

(32a) 



— {PaP - pSap) +PaP 
Tr 



,,^dT{t) 1 
-L{t)—^ qa+qa. 

dra Tr 



(32b) 



To obtain the hydrodynamic balance equations, we 
substitute Pa/? and qa into equations l|5b(l and (|3c|l . which 
are the conservations laws. First we calculate some of the 
essential expressions that will be needed later: 



dpap ^ ^_J^2^ +i_^(v.M) 
drp drp \ " 3 drp 



dfif Idu 

-2 Aap - 

rp \ 3 dry 

^ _9 / dpap _ dp_ 
^ dt\ drp drp 



'5, 



ap 



(38) 



where 



Pap 



dt 



dua _^dup 2 duy 
drp dra 3 dr^ 



5ap^ dt' , (33a) 



PaP 



dUa 

dra 



r.^+pj(V.«)-2MA.,^ + 

2 2 dpap dUa , . 

+ -MV.«) (39) 



dt dra 
are the last terms in Eqs. ((32(1 . 

dM{t,t') m f fdu. 



dt 
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dug 

dt ^/' + ^^")^ 



(34a) 



^^-XV^T-VXWT-Tr-^^. (40) 
dra dt dra 

Taking into account l(5b(l and ((38(1 we obtain: 



dt 



( dua 



dUa 

drp 



( dUa 



dUa 

drp 



dra \ 3 



^2 ^P 



(41) 
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where 



Z = 2 



9/i 

drr. 



A 1 ^""7 r 
Aq,3 - 0(1/3 

3 or. 



From Eqs. |ScI), and (gUl one has 



/9T 



dt 



- up 



dT 
drp 



dT 
'dt 



up 



dT 
drp 



AV^T + VAVT - 

dp \ /„ X „ , dua 2 ,„ ,o 

^r-^+pj (V-IX) + 2^mKp— - -M(V-lt)2 



9f i9r^ 



d_ 

^dt 



dp 
"dt 



dUg 



gpaff dUg 

dt drp 



In Eqs. H41(l and 142|l terms of the first order are fi, A, 
u, du/dxi, dp/dxi and dT/dxi. We keep only the terms 
of the first order and neglect all terms which contain the 
derivatives over and A. Then in Eq. H41|) one can 
neglect Z, and, in the right-hand part of Eq. (|42l) . the 
second, fourth, fifth term, and, in the square brackets, 
the second and third term. For the sixth term, and the 
last term in the square brackets on the right-hand side of 
Eq. (|42l) one can apply the following approximation: 



dpap dug 
dt drp 

dkap 



dp 
'dt 



5a 



'dt 



A, 



1 



aP 



-2p. 



dt 



d^ 



-rrg^iP^P 



■ pSap) 



dUa 

drp 



dp 
dt 



(V-«) 



d'^Pap dUg 
^ df^ drp 



d f dp , 

»|(v.«,, 



d'^Pap dUg 
dt"^ drp 



(43) 



where we use the same above mentioned rules of the first 
order approximation. Then we obtain the hyperbolic hy- 
drodynamic balance equations of the first order 



| + V(p.) = 0, 
the equation of continuity. 



(44) 



p(^ + (V.«)«) 



/iV^M (45) 



the hyperbolic Navier-Stokes equation. 



9 , 

^^■dt^' 



mX ,„ 1 , , / d 



-V^T--(V.«)^..- + lj(pT) (46) 

— the hyperbolic heat conduction equation. Here, c„ = 
3/2 and p — nkT. 

If ti = then Eq. (|^ takes the form: 



d^T dT mX 



dt^ dt Cykp 



(47) 



(42) 



which is the hyperbolic heat conduction equation in static 
medium. 

Equations 144|l - 147() are transformed to the classical 
hydrodynamic balance equations if = 0. 



V. HYPERBOLIC HEAT FLUX AND 
HYPERBOLIC HEAT CONDUCTION EQUATION 

In this section we consider the heat propagation in 
static medium, i. e. when the mass velocity equal zero 
u = 0, based on the results obtained for the heat fiux and 
the heat conduction equation in the previous section. 

One can see that Eq. Ij36b|) is the hyperbolic heat flux 
Eq. where the hyperbolic relaxation time Tg is the 
Maxwellian relaxation time r„: 



-X 



dT 
drg 



dqg 

dt ' 



(48) 



where the thermal conductivity A is defined by Eq. Ij30b|l . 
Thus, we show that the Cattaneo equation is the 
result of the relaxation time approximation. 

In the limit — > we obtain the classic form for the 
heat flux 



lim qg 



-A 



dT 
drg ' 



(49) 



which corresponds to the case when t ^ Tr (cp. Eq. 
(I29HH . 

Let us write the hyperbolic heat conduction equation 
or telegraph equation l|T7|) : 



where 



d^T dT 



mA 

Cvkp 



(50) 



is the thermal diffusivity. Taking into account that the 
specific heat is C = Cyk/m, the thermal diffusivity is 
written as 



A 



(51) 
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which is in agreement with the classic definition 0, 125I 

Denoting = (5/6)u^ and substituting A from l)30b|l 
we obtain the expression for the a: 



a = TrW 



From Eq. (|50|l one can see, that for 
the parabohc heat conduction equation 



dT 
'dt 



drl ' 



(52) 
we have 



(53) 



For Tr ^ 00 the expression ()5U|) becomes a wave equa- 
tion 



drl ' 



(54) 



whose solutions are known in the literature as the second 
sound, which was observed in solids at low temperatures 

The relaxation time is the time of establishing of the 
local equilibrium Maxwell distribution (local thermody- 
namical equilibrium) by the definition from Eq. Q. The 
exact value for the depends on the type of the molec- 
ular interaction, i.e. on the collision integral (cp. Eqs. 
(O and O). It is known that the relaxation time Tr and 
the mean free time tq have the same order. However, let 
us estimate it on the other hand. 

As follows from structure of Eq. 15UI) the heat front is 
propagated with the speed Vf 4, 10, 35, 36] : 



(55) 



i.e. practically with the average molecular velocity. If we 
suppose that the distance of the propagation of the heat 
front for the time is not less than the mean free path 

TrVf > TqV, (56) 

the estimation for the relaxation time has the form 



(57) 



One can see that the expression H57|l is not in distinction 
with the known fact that the equilibrium Maxwell distri- 
bution is established over 3-4 molecular collisions l25Ll2^ . 
In this way, the equations (|45|) - (|47|l are suitable for the 
description of the processes at time intervals t which have 
order of the mean free time tq for molecules. The limit 
To ^ corresponds to a process when t ^ tq, i. e. when 
the processes described by the classical Navier-Stokes and 
heat conduction equations were already formed. 

In conclusion of this section we estimate the relax- 
ation times for Nitrogen and phonons in Aliminium. The 
expression for the Maxwellian relaxation time of mono- 
atomic gas is very simple (see Eq. I30a() : 



p 



(58) 



where p is the static pressure. The relaxation time for 
phonons one can obtain in the approximation of phonon 
gas ^ : 



X = ^CpulTr 



(59) 



where Us is the sound speed. Substituting /i=4.15 Poise 
for Nitrogen at 1000 K and atmospheric pressure [ssj . 
and A=209.3 W/m K, C=880 J/kg, p=2.7 g/cm^ and 
Us=5.8 km/s for Aluminium [33 we obtain for Tr the 
values 4.096x10-1° s and 1.024x10"" s for Nitrogen 
and Aluminium, respectively. As one can see from the 
comments to Eq. these values for the Maxwellian 
relaxation times are in a good agreement with the relax- 
ation times for the hyperbolic heat conduction given by 
Lykov [3. 



VI. CONCLUSION 

In the present work, we have applied the integral form 
for the Boltzmann equation and the relaxation time ap- 
proximation for the collision integral. In the first order 
iteration we have derived the hyperbolic-like forms for 
the heat transfer and for the Navier-Stokes equations. 
For the static medium we have obtained the Cattaneo 
heat flux and the well known telegraph equation for the 
heat conduction. It was shown that the relaxation time 
for the hyperbolic equations is the Maxwellian relaxation 
time. The obtained equations can be used for the descrip- 
tion of the fast heat-mass transfer processes of the order 
of the mean free time. 

Additionally, for the heat propagation in static 
medium, we estimate the hyperbolic and Maxwellian re- 
laxation times for the heat conduction in Nitrogen and 
Aluminium. The obtained values are found to be in a 
good agreement with experiment. 



Appendix A: CALCULATION OF M{t, t) AND L{t, t) 

1. 

Let us write the expression (|28a|) for M{t) — M{t,t): 



M{t) = — / mcacp 



( 1 2 e 



d?c, (Al) 



M{t) = 



lORT 



(0) 



(CaC/3 ) (Cq, C0 ) - - (Cq C/5 ) (5q^ ) 



(A2) 



For the expression in brackets we have: 

1 2 , ^ ^ ^ 
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Then M takes the form: 

[3=1 a=l 

and it can be separated into two parts: 

M(i) = 3J„„ + 6J„^, (A5) 
where Jaa and Jq/j have simple expressions: 

j^^ ^ [ fWcU'c = -nkT, (A6) 

10i?r3y ^ c„c^d c- ^^nfci. 
Finally, Af has the form: 

M(i) = nkT (A8) 

2. 

Let us write now the expression (|28b|l for L{t) = 
L{t,t): 
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